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. Abstract. On the basis of the competing cubic-quintic nonlincarity model, stability (insta- 

CN| ' bility) of continuous waves in nonlocal random non-Kcrr nonlinear media is studied analy- 

tically and numerically. Fluctuating media parameters are modeled by the Gaussian white 
^ ' noise. It is shown that for different response functions of a medium nonlocality suppresses, 

. as a rule, both the growth rate peak and bandwidth of instability caused by random pa- 

rameters. At the same time, for a special form of the response functions there can be an 
"anomalous" subjection of nonlocality to the instability development which leads to further 
increase of the growth rate. Along with the second-order moments of the modulational 
amplitude, higher-order moments are taken into account. 
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1 Introduction 

> 

Modulational instability (MI) is a universal phenomenon arising in many nonlinear systems. It is 
closely related to formation of localized wave structures - solitons and considered as a precursor 
of this formation. In optics, in case of Kerr-like nonlinear media MI is responsible for formation 
of bright solitons and its absence is necessary for dark solitons formation pQ . MI consisting in the 



OO 

' exponential growth of small perturbations in the amplitude or the phase of initial optical waves 

leads to the breakup of a beam or a quasi-coherent pulse into filaments that can evolve into 
trains of solitons [2]. MI finds its application in plasma [3], hydrodynamics [4], fluids [5], atomic 
Bose-Einstein condensates [H [7] . The realm of applications also includes optical communications 
systems, all-optical logical devices [8], and the generation of ultra-high repetition-rate trains of 
soliton-like pulses [9]. 

In the setting of nonlinear optics it was shown that MI is absent for defocusing Kerr-type 
media, the medium being described by deterministic parameters. In case of focusing determinis- 
tic Kerr media MI presents as the long- wave instability with a finite bandwidth [10]. In frames 
of a more realistic approach that implies random behavior of the characteristic parameters of 
a medium, that is, they are considered to fluctuate randomly around their mean values, the 
region of MI was shown to be extended by stochastic inhomogeneities of a Kerr medium over 
the whole spectrum of modulation wave numbers. This is true for both focusing and defocusing 
regimes 

The next step toward more realistic description implies taking into account nonlocality of 
a medium. Moreover, some media manifestly exhibit nonlocal properties [12] and should be 
described by appropriate nonlocal models. As a rule, this property is a result of a corresponding 
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transport process such as heat conduction in thermal nonlinear media |13j . diffusion of atoms in 
a gas |14j , many-body interaction in Bose-Einstein condensates [TJ]. Nonlocality, as compared 
to local cases, can crucially change the behavior of proceeding in such a medium processes. For 
example, it can support the propagation of self-focused beams without collapsing [16l[T7]. Dark 
solitons interaction [18] is greatly influenced by nonlocal properties of a medium. The investi- 
gation of MI in deterministic nonlocal Kerr-type media revealed its absence for the defocusing 
case for small and moderate values of the product "modulation amplitude x nonlocality pa- 
rameter" [19] . The overview of the modulational instability and beam propagation in nonlocal 
nonlinear media can be found in |20j. 

The incorporation of nonlocality into a model under study is generally accomplished via the 
nonlinear refractive index of a material, whose dependence on the intensity is determined by the 
chosen model of nonlinearity. Usually, experiments show the deviation from the linear (Kerr) 
dependence of the refractive index on the incident intensity for large intensities. This occurs 
for a range of materials: semiconductor-doped glasses [21 j . semiconductor waveguides [22] . and 
nonlinear polymers [23]. Saturable and cubic-quintic nonlinearities are the most typical for 
nonlinear systems. 

The objective of this paper is to investigate MI in nonlocal medium with stochastic param- 
eters and competing cubic-quintic nonlinearity. The fact that nonlocality spreads out localized 
excitations allows to anticipate the decrease of both the growth rate peaks and bandwidths of 
instability. Such a situation is valid for the case of stochastic media with the sign-definite Fourier 
images of the response functions. However, for nonlocal media with sign-indefinite Fourier image 
of the response function MI gain can exceed that of a local stochastic medium. This behavior 
can be considered as "anomalous" . 

In the present paper the investigation of MI of continuous waves in nonlocal stochastic media 
is based on a generalized nonlocal nonlinear Schrodinger equation with random coefficients and 
competing nonlinearities. For the illustration of the obtained results we use the white noise 
model for parameter fluctuations and response functions of several types. 



2 Model 

In this paper we consider a medium with nonlinearity that is a nonlocal function of the incident 
field. This nonlinearity induced by a wave with the intensity I(x, z) can be presented in general 
form 

/oo poo 
dx'R(x-x')I a (x',z) + s(z) i dx'L(x- x')I 2a (x',z), (1) 
-oo J — oo 

where x is the transverse coordinate, a is a positive constant, z can be considered as time or 
the spatial propagation coordinate depending on the context, nonlinearity coefficients g(z) and 
s(z) are considered as stochastic functions which fluctuate around their mean values go and sq: 

g{z) = go{l+m g (z)), s(z) = s (l + m s (z)), (2) 

Here m g and m s are independent zero-mean random processes of the Gaussian white-noise type, 

(m g ) = (m s ) = 0, (m g (z)m g (z')) = 2a g 5(z - z'), (m s (z)m s (z')) = 2a 2 s b(z - z'). 

The angle brackets designate the expectation with respect to the distribution of the correspon- 
ding processes. In the imposed phenomenological model ([1]) the field-intensity dependent change 
of the refractive index is characterized by two normalized symmetric response functions R(x) 
and L(x), J^° oQ dxR(x) = J^° oo dxL(x) = 1. In case when R(x) = L(x) = 5(x) (local limit) and 
constant nonlinearity coefficients we reproduce the well-known general power-law dependence 
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on the incident intensity for local models with competing nonlinearities pQ. At the same time, 
for equal response functions one can consider the expression as the second-order expansion of 
the model with saturable nonlinearity, provided a = 1 and s = —g/I s = —rioc/I 2 , where I s is 
the saturation intensity and corresponds to the maximum change in the refractive index. 
The set of the parameters a = 1 and s(z) = is associated with nonlocal Kerr- type media with 
random parameters [24j . 

Thus the propagation of plane waves along the z axis in a stochastic medium with nonlocal 
competing nonlinearities is governed by the generalized nonlinear Schrodinger equation 

iu z + ^d(z)u xx + uAn(\u\ 2 ) = 0, (3) 

where u(x, z) is the complex envelope amplitude, A?i(|u| 2 ) is defined by (pQ) and standard dimen- 
sionless variables are used. The diffraction coefficient d(z), just as g{z) and s(z), is a stochastic 
function with the same properties: 

(m d ) = 0, (m d (z)m d (z')} = 2a 2 d 5{z - z'). 

Equation ([3]) represents the generalization of the Kerr-type nonlocal model considered in [23] for 
the case of higher-order nonlinearities. The deterministic variant was comprehensively studied 
in [25]. 

In further considerations we adopt that do > and the condition go • so < stands for the 
competition of nonlinearities. 

As a solution equation © admits a homogeneous plane wave 



uq = A exp 



iA 2a [ Z dz' (g(z')+A 2a s(z')) 
Jo 



(4) 



where A is a real amplitude. Proceeding to a linear stability analysis of the solution d3|) we 
assume that 



u(x, z) = (A + v(x, z)) exp 



iA 2a / dz' (g(z') + A 2a s{z')) 



ii 



(5) 



is a perturbed solution of equation ([3]). Here v(x,z) is a small complex modulation. The 
substitution of equation ([5]) into equation ([3]) and the linearization around the plane wave dH) 
provide a linear equation for v(x,z): 

iv z + ^d(z)v xx + 2ag(z)A 2a J dx R(x - x')Re v(x', z) 

+ 4„ 5W ^/d^-,OR-(x',z) = 0. (0) 

Decomposing v into real and imaginary parts, v = v r (x,z) + ivi(x,z), and performing the 
Fourier transforms 

2 poo 2 roo 

p(k,z) = — dxv r (x,z)e lkx , a(k,z) = — dx Vi(x, z)e lkx , 

i poo i roc 

R(k,z) = — dxR(x,z)e ikx , L(k,z) = — dx L(x, z)e ikx , 
2k J-oo 2ir 



we convert equation ([6]) to a system of linear equations for p and a: 



d ( P \ _ ( \d(z)k 2 \ ( p 



dz\ a J y -\d{z)k 2 + 2ag{z)A 2a R + 4as{z)A Aa L J V a J ' ^ 

Equation ([7]) is sufficient for the investigation of MI for a deterministic medium with the 
parameters do, go, and so- in contrast, for stochastic systems one should consider the second 
moments (p 2 ), {per) and (a 2 ) providing the minimal nontrivial information about MI induced 
by random fluctuations. 
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3 The second-order moment MI gain 

3.1 Theoretical approach 

We consider the vector of the second moments 

X^ = ((p 2 ),(pa),(a 2 )f. (8) 
The Furutsu— Novikov formula 



( Ui (r)F[v])=J2j dr/ ^( r - r ')(^)' (9) 

where u is an arbitrary zero-mean Gaussian field with the covariance function Bij(r), F[u] is 
a functional, allows to express products like (m^p 2 ) in terms of components of X^ 2 \ As a result, 
is subject to the evolution equation (d/dz) X^ = M^X^ 2 \ where M^> is the 3x3 matrix 
of the form 



= \ C - \B -A \B , (10) 





where 

A = <r 2 d d 2 k A , B = d k 2 , C = 2ag A 2a R + 4as A Aa L, 
D = 6Aa 2 a 2 s 2 A 8a L 2 + 16a 2 a 2 g g 2 A ia R 2 . 

Instability occurs for positive real parts of the eigenvalues of M^ 2 ' . MI gain G2(k) is determined 
by the largest positive value. We will illustrate the results by means of the exponential response 
function with the sign-definite Fourier image 

*.(«)- i«p(-M), *.(*)- (id 

and as an example of the sign-indefinite Fourier transform response function we will take the 
rectangular response function 

/ \ 7TT for |x| < A, - sin(Afc) . , 

R r (x) = l 2A Rr{k) = ^_±, ( 12 ) 

^ for \x\ > A, 

where A is the nonlocality parameter. Below we will study in detail the case of competing cubic- 
quintic nonlinearity, that is, a = 1, go > and sq < and the response functions are assumed 
to be equal. For analysis purposes it is convenient to introduce a quantity 

x - 



2\s \A 2 



with the help of what we will discriminate between the focusing and the defocusing regimes. 
The inequality x > 1 corresponds to the predominance of the cubic nonlinearity term over the 
defocusing one, x < 1 determines the defocusing regime when the quintic term prevails. 
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Figure 1. Plots of the MI gain G*2(fc) for: (a) local stochastic media (solid lines) and nonlocal stochastic 
media with the exponential response function (dashed lines) for different values of the parameter \: the 
pair of red curves corresponds to \ — 2, blue ones - to x — lj green ones - to % = 0.8, here A = 1; (b): 
local stochastic media (blue line) and nonlocal stochastic media with the exponential response function 
for x = 0.8: A 2 = 1 (red line), A 2 = 18 (green line); (c): local deterministic media with competing 
cubic-quintic nonlinearity (yellow line), local stochastic media (blue line), and nonlocal stochastic media 
with the exponential response function for x = 1-4: A 2 = 1 (red line), A 2 = 18 (green line). On all plots: 
d = 2, g A 2 = 1, a 2 = a 2 = a 2 = 0.1. 
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Figure 2. Plots of the MI gain G2(k) for local stochastic media (blue lines) and nonlocal stochastic 
media with the rectangular response function: (a) x = 0.8, (b) x = 1-4. On both plots A = 2 (red lines), 
A = 6 (green lines), A = 25 (yellow lines), d = 2, g A 2 = 1, a 2 = a 2 = a 2 = 0.1. 



3.2 Numerical results 

In Fig. QJa) we demonstrate the influence of the transition from the focusing regime (x > 1) 
to the defocusing one (% < 1) on the MI bandwidth and the MI gain peak for media with 
cubic-quintic nonlinearity as the intensity of the incident field is increased. This transition is 
accompanied by the pronounced enhancement of the MI bandwidth. At the same time the MI 
gain peak is decreased. For the defocusing regime the matrix (|10p possesses one real eigenvalue Ai 
that is positive for all k and two complex conjugate ones A2 and A3 with negative real parts 
in the case of the exponential response function. The growth rate peak of G%(k) = Ai and MI 
bandwidth are suppressed by nonlocality. The growth of the nonlocality parameter A results in 
the amplification of the suppression effect (Fig.rjjb)). The situation for the rectangular response 
function f)12|) is different. For sufficiently high nonlocality MI gain maximum for a given wave 
number k can exceed the corresponding value of G2 for a local random medium with cubic-quintic 
nonlinearity (left panel of Fig. [2j). At the same time, the MI bandwidth becomes strictly finite 
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in this limit A local deterministic medium with the competing cubic-quintic nonlocality in the 
focusing regime (x > 1) (Fig. Die)) produces the long- wave instability with a finite bandwidth. 
The bandwidth is extended by stochasticity of medium parameters to the whole spectrum of 
modulation wave numbers. Corresponding calculation of eigenvalues of the matrix M2 (|10p for 
the case of x > 1 demonstrates the suppression the MI gain and bandwidth for media with both 
sigh-definite and sign- indefinite response functions (Fig. Efb)). The positions of the MI gains 
shift toward smaller wave numbers k under nonlocality growth, producing finite bandwidth. 

4 Higher-order moments 

4.1 Theoretical approach 

Evidently, in frames of the second-order moments analysis (JH) it is impossible to characterize 
in full detail the dynamics of a stochastic system. For example, featured points of MI gain 
(maximum position, bandwidth, etc.) calculated from the second moments could fluctuate 
when accounting for higher-order moments. More detailed information about MI is provided by 
the higher-order moments 

X^) = {{p^-i) a i)Y i = 0,...,2n. (13) 
As before, applying the Furutsu-Novikov formula, we obtain a matrix M^ 2n ^ in the form 

= d k 2 A^ + (2ag A 2a R + 4:asoA 4a L - ±d k 2 )B^ 

+ d 2 k 4 a 2 d C^ + {6Aa 2 s 2 A 8a La 2 + 16a 2 g 2 A 4a R 2 a 2 )D^ (14) 

with the following non-zero entries of the matrices A^ 2n \ B^ 2n \ and L>( 2n ) 

^=(»-Q(»-^)> cgi-flW-i^-D, j u ..... 

Then among the roots of the characteristic polynomial det |M^ 2n ) — XI \ the maximal real part 
will determine nG2 n (k)- The fact that the characteristic polynomial is of the odd degree and all 
the matrix elements of M^ 2n ^ are real ensures at least one real eigenvalue of M^ 2n \ the others 
being mutually complex conjugate. Below the analysis of the 4-th and 6-th moments is given. 

4.2 Numerical results 

Fig. [3] presents the results of calculating MI gains G2, G4 and Gq for the exponential response 
function. In case of the focusing regime as well for the defocusing one nonlocality suppresses 
the higher-order moments. For the focusing regime the maxima of MI gains of higher orders are 
shifted to higher wave numbers k. The opposite shift direction occurs for the defocusing regime. 
As compared to the case of Kerr- type media |24j , this apparently exhibits the stochastic origin of 
the MI gains. The similar situation takes place for media with the rectangular response function 
(Fig. 2£ a)). Fig. IHb) provides the dependence of the maximum of the main MI gain peak on 
the nonlocality parameter A for the case when the quintic nonlinearity term dominates (\ < !)• 
Such a dependence reflects the possibility for MI gain maximum to exceed the corresponding 
value of MI gain for a local random medium in a narrow domain of wave numbers k (see 
Fig. [2]^a)). At the same time Fig. H^b) depicts the closing in of MI gains of different orders as 
the nonlocality is increased, that is, high nonlocality smooths the fluctuations of the modulation 
amplitude growth. This is valid irrespectively of the response properties of the media regarded 
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Figure 3. Plots of the MI gains Ge(k) (solid lines), G4 (dashed lines), and G2 (dash-dotted lines) for 
local stochastic media (red curves on both plots) and nonlocal stochastic media with the exponential 
response function; (a): x = 2, A = 2 (yellow curves), A = 4 (green curves); (b): \ = 0.9, A = v2 (yellow 
curves), A = 4 (green curves). On both plots do = 2, goA 2 = 1, u\ = a 2 = erf = 0.1. 




Figure 4. (a) Plots of the MI gains Ge(k) (solid lines), G4 (dashed lines), and G2 (dash-dotted lines) for 
local stochastic media (blue curves) and nonlocal stochastic media with the rectangular response function 
(red curves); here A = 1. (b) Plots of the maximum of the main peak of the MI gains G§(k) (green line), 
G4 (red line), and G2 (yellow line) versus nonlocality parameter A. On both plots x = 0-8, do = 2, 
g A 2 =l >a 2 d = a 2 g = a* = 0.1. 



and regimes of cubic-quintic nonlinearity. In Fig. [5] we demonstrate MI gains G2, G4 and Gq for 
the exponential response function for the regimes with diverse signs of nonlinearity coefficients 
than considered above and compare them with the competing cubic-quintic nonlinearity regime. 
For focusing cubic and quintic nonlinearities (go > 0, sq > 0) one observes the same characteristic 
pattern consisting in the suppression MI gains, at the same time as compared to competing 
nonlinearity MI gains are less suppressed (Fig. GD(a)). The opposite situation takes place for the 
case of defocusing cubic and focusing quintic nonlinearities (go < 0, so > 0). The MI gain peaks 
are suppressed more than their competing cubic-quintic nonlinearity counterparts (Fig. 0(b)). 
In both cases one can see the transition from the initial regime (initially defocusing regime on 
Fig. GD[a) and focusing one on Fig. EJ^b)) as one changes the chosen set of signs of nonlinearity 
coefficients. 
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Figure 5. Plots of the MI gains G§(k) (solid lines), G4 (dashed lines), and G2 (dash-dotted lines) for 
nonlocal stochastic media with the exponential response function for competing cubic-quintic nonlincarity 
with g > and s < (red curves on both plots) and: (a) go > and so > (green curves), here \ = 0-9; 
(b) go < and so > (green curves), here \ = 2- On both plots do = 2, g$A 2 = 1, u\ = er^ = of = 0.1, 
A = 1. 

5 Conclusion 

In this paper for the phenomenological model supporting two regimes of the wave propagation, 
applying the linear stability analysis we have investigated the MI of a homogeneous wave in 
a nonlocal non-Kerr medium with cubic-quintic nonlinearity and random parameters. Adopting 
a white-noise model for parameter fluctuations, we have obtained the equations determining the 
dependence of the MI gain on the modulation wave number. Nonlocality was proved to suppress 
considerably the stochasticity-induced MI growth rate for media with the sign-definite Fourier 
images of the response functions. At the same time, it was shown that for the nonlocal media 
with the sign-indefinite Fourier images of the response functions MI gain can exceed that for 
local media for some wave number intervals. Some attention was paid to the cases of focusing 
cubic and quintic nonlinearities and defocusing cubic and focusing quintic nonlinearities. 
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